n_distinct(result$MA.id) # number of MAs
## [1] 23492
nrow(result) # 9 model result per each MA
## [1] 211428
table((result %>% distinct(MA.id, .keep_all = TRUE))$event)
##
## common infrequent rare very rare
## 22074 1154 101 163
proportions(table((result %>% distinct(MA.id, .keep_all = TRUE))$event)) * 100
##
## common infrequent rare very rare
## 93.9639026 4.9123106 0.4299336 0.6938532
There are three types of model that have \(\tau\)
ivre,
htevaguelogit and htelogit_tau50 (log
ratios)ablogit and
ablogit_tau50 (log odds)htebeta
(prob)cor(
result_bayesian[result_bayesian$modelname == "htelogit_tau50",]$tau,
result_bayesian[result_bayesian$modelname == "ablogit_tau50",]$taua
)
## [1] 0.2477361
cor(
result_bayesian[result_bayesian$modelname == "htelogit_tau50",]$tau,
result_bayesian[result_bayesian$modelname == "ablogit_tau50",]$tauc
)
## [1] 0.3046824
cor(
result_bayesian[result_bayesian$modelname == "htebeta",]$taua,
result_bayesian[result_bayesian$modelname == "ablogit_tau50",]$taua
)
## [1] -0.08585503
cor(
result_bayesian[result_bayesian$modelname == "htebeta",]$tauc,
result_bayesian[result_bayesian$modelname == "ablogit_tau50",]$tauc
)
## [1] -0.1115232
| modelname | Sig_sum | percentage_total |
|---|---|---|
| ivfe | 10403 | 44.283160 |
| ivre | 9132 | 38.872808 |
| ctevaguelogit | 11553 | 49.178444 |
| ctebeta | 11402 | 48.535672 |
| htevaguelogit | 7308 | 31.108462 |
| htelogit_tau50 | 6047 | 25.740678 |
| htebeta | 4363 | 18.572280 |
| ablogit | 8003 | 34.066916 |
| ablogit_tau50 | 597 | 2.541291 |
table(result_sigsum$Sig_sum)
##
## 0 1 2 3 4 5 6 7 8 9
## 10947 1075 706 1407 1587 922 915 2312 3121 500
round(proportions(table(result_sigsum$Sig_sum))*100,3)
##
## 0 1 2 3 4 5 6 7 8 9
## 46.599 4.576 3.005 5.989 6.755 3.925 3.895 9.842 13.285 2.128
table(result_sigsum$event, result_sigsum$Sig_sum)
##
## 0 1 2 3 4 5 6 7 8 9
## common 9948 977 632 1340 1524 902 883 2261 3109 498
## infrequent 799 74 58 59 61 18 29 42 12 2
## rare 65 8 7 5 2 2 3 9 0 0
## very rare 135 16 9 3 0 0 0 0 0 0
round(proportions(table(result_sigsum$event, result_sigsum$Sig_sum), margin = 1)*100,2)
##
## 0 1 2 3 4 5 6 7 8 9
## common 45.07 4.43 2.86 6.07 6.90 4.09 4.00 10.24 14.08 2.26
## infrequent 69.24 6.41 5.03 5.11 5.29 1.56 2.51 3.64 1.04 0.17
## rare 64.36 7.92 6.93 4.95 1.98 1.98 2.97 8.91 0.00 0.00
## very rare 82.82 9.82 5.52 1.84 0.00 0.00 0.00 0.00 0.00 0.00
For best model comparison among seven Bayesian models (exclude ivfe and ivre)